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. _ ' Abstract. 

^nq ■ In this paper we set bounds on the radiation content of the Universe and neutrino 

properties by using the WMAP-5 year CMB measurements complemented with most 
' (-('' of the existing CMB and LSS data (WMAP5+A11), imposing also self-consistent BBN 

^ I' constraints on the primordial helium abundance. 

^ I We consider lepton asymmetric cosmological models parametrized by the neutrino 

\-^ ' degeneracy parameter ^i, and the variation of the relativistic degrees of freedom, 

C/3 , AN°f^f , due to possible other physical processes occurred between BBN and structure 

formation epochs. 

We get a mean value of the effective number of relativistic neutrino species of 
Neff = 3.256^Qg4^ (68% CL), bringing an important improvement over the similar 
result obtained from WMAP5-fBA0-|-SN-|-HST data [57]. 
C"*~- . We also find a strong correlation between Vl„ih? and z^q, showing that we observe N^ff 

mainly via the effect of Zeq^ rather than via neutrino anisotropic stress as claimed by 
the WMAP team [57]. 

WMAP5-I-A11 data provides a strong bound on helium mass fraction of 1^ = 
QQ ' 0.2486 ± 0.0085 (68% CL), that rivals the bound on Yp obtained from the conservative 

^^ ■ analysis of the present data on helium abundance. 

For neutrino degeneracy parameter we find a bound of —0.216 < ^j/ < 0.226 
(68%CL), that represent an important improvement over the similar result obtained 
by using the WMAP 3-year data. 
C^ \ The inclusion in the analysis of LSS data reduces the upper limit of the neutrino 

mass to m^ < 0.419 eV (95% CL) with respect to the values obtained from the analysis 
from WMAP5-only data [56] and WMAP5-hBA0-t-SN-fHST data [57] . 

We forecast that the CMB temperature and polarization measurements observed 
with high angular resolutions and sensitivities by the future Planck satellite will 
reduces the errors on ^i, and Yp down to (j{£^v) ~ 0.089 (68% CL) and aiYp) = 
0.013 (68% CL) respectively, values fully consistent with the BBN bounds on these 
parameters. 

This work has been done on behalf of Planck-LFI activities. 



PACS numbers: CMBR theory, dark matter, cosmological neutrinos, big bang 
nucleosynthesis 
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1. Introduction 

The radiation budget of the Universe relies on a strong theoretical prejudice: apart from 
the Cosmic Microwave Background (CMB) photons, the relativistic background would 
consist of neutrinos and of possible contributions from other relativistic relicts. The main 
constraints on the radiation energy density come either from the very early Universe, 
where the radiation was the dominant source of energy, or from the observation of 
cosmological perturbations which carry the information about the time equality between 
matter and radiation. 

In particular, the primordial light element abundance predictions in the standard 
theory of the Big Bang Nucleosynthesis (BBN) [U El [3l H] depend on the baryon-to- 
photon ratio, rjB, and the radiation energy density at the BBN epoch (energy density 
of order MeV^), usually parametrized by the effective number of relativistic neutrino 
species, Ngff. 

Meanwhile, the number of active neutrino flavors have been fixed by Z^ boson decay 
width to Ni, = 2.944 ± 0.012 |4] and the combined study of the incomplete neutrino 
decoupling and the QED corrections indicate that the number of relativistic neutrino 
species is N^jf = 3.046 [5]. Any departure of A^e// from this last value would be due to 
non-standard neutrino features or to the contribution of other relativistic relics. 
The solar and atmospheric neutrino oscillation experiments [6|, [7] indicate the existence 
of non-zero neutrino masses in eV range. 

There are also indications of neutrino oscillations with larger mass-squared difference, 
coming from the short base-line oscillation experiments [H [9], that can be explained 
by adding one or two sterile neutrinos with eV-scale mass to the standard scheme with 
three active neutrino flavors (see Ref . jlOj for a recent analysis) . Such results have impact 
on cosmology because sterile neutrinos can contribute to the number of relativistic 
degrees of freedom at the Big Bang Nucleosynthesis [11]. These models are subject to 
strong bounds on the sum of active neutrino masses from the combination of various 
cosmological data sets [l2l[T3], ruling out a thermalized sterile neutrino component with 
eV mass [TitfTS]. 

However, there is the possibility to accommodate the cosmological observations with 
data from short base-line neutrino oscillation experiments by postulating the existence of 
a sterile neutrino with the mass of few keV having a phase-space distribution significantly 
suppressed relative to the thermal distribution. 

For both, non-resonant zero lepton number production and enhanced resonant 
production with initial cosmological lepton number, the keV mass sterile neutrino 
produced via small mixing angle oscillation conversion of thermal active neutrinos 
[16] provides a valuable Dark Matter (DM) candidate, alleviating the accumulating 
contradiction between the ACDM model predictions on small scales and observations 
by smearing out the small scale structure [171 [THl UHl [20] . 
A recent analysis of X-ray and Lyman-a data indicate that keV sterile neutrinos can be 
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considered valuable DM candidates only if they are produced via resonant oscillations 
with non-zero lepton number or via other (non-oscillatory) production mechanisms |21j . 
On the other hand, the possible existence of new particles such as axions and gravitons, 
the time variation of the physical constants and other non-standard scenarios (see e.g. 
[22] and references therein) could contribute to the radiation energy density at BBN 
epoch. 

At the same time, more phenomenological extensions to the standard neutrino 
sector have been studied, the most natural being the consideration of the leptonic 
asymmetry [23l [2ll [25] , parametrized by the neutrino degeneracy parameter ^j, = fiu/T^o 
[fiu is the neutrino chemical potential and T^^^ is the present temperature of the neutrino 
background, T,JT,^^ = (4/ll)i/3]. 

Although the standard model of particle physics predicts the value of leptonic 
asymmetry of the same order as the value of the baryonic asymmetry, B ~ 10~^°, there 
are many particle physics scenarios in which a leptonic asymmetry much larger can be 
generated [261 [27] . One of the cosmological implications of a larger leptonic asymmetry 
is the possibility to generate small baryonic asymmetry of the Universe through the 
non-perturbative (sphaleron) processes [281 [291 [30] . Therefore, distinguishing between 
a vanishing and non-vanishing C,u at the BBN epoch is a crucial test of the standard 
assumption that sphaleron effects equilibrate the cosmic lepton and baryon asymmetries. 
The measured neutrino mixing parameters implies that neutrinos reach the chemical 
equilibrium before BBN [311 l32l [33] , so that all neutrino flavors are characterized by 
the same degeneracy parameter, ^,y, at this epoch. The most important impact of the 
leptonic asymmetry on BBN is the shift of the beta equilibrium between protons and 
neutrons and the increase of the radiation energy density parametrized by: 



AAT.^Xe.) = 3 



30 /e.y 15 /e. 



(1) 



The BBN constraints on Ng// have been recently reanalyzed by comparing the 
theoretical predictions and experimental data on the primordial abundances of light 
elements, by using the baryon abundance derived from the WMAP 3-year (WMAP3) 
CMB temperature and polarization measurements [3ll[35l[36]: t]b = 6.14 x 10^^°(1.00± 
0.04). In particular, the ^He abundance, Yp, is quite sensitive to the value of A^e//- In 
the analysis of Ref. [37], the conservative error of helium abundance, Yp = 0.249 ±0.009 
[38] , yielded to A'e// = 3.1;^};2 (95% CL) in good agreement with the standard value, 
but still leaving some room for non-standard values, while more stringent error bars of 
helium abundance , Yp = 0.2516 ±0.0011 [39], leaded to N^ff = 3.32+°:^3 (95% CL) ^. 

The stronger constraints on the degeneracy parameter obtained from BBN [H] 
gives —0.04 < ^ < 0.07 (69% CL), adopting the conservative error analysis of Yp of Ref. 
[38] and ^ = 0.024 ± 0.0092 (68% CL), adopting the more stringent error bars of Yp of 
Ref. 



The CMB anisotropics and LSS matter density fluctuations power spectra carry the 
signature of the energy density of the Universe at the time of matter-radiation equality 
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(energy density of order eV^), making possible the measurement of N^ff through its 
effects on the growth of cosmological perturbations. 

The number of relativistic neutrino species influences the CMB power spectrum by 
changing the time of matter-radiation equality that enhances the integrated Sachs- Wolfe 
effect, leading to a higher first acoustic Doppler peak amplitude. Also, the temperature 
anisotropy of the neutrino background (the anisotropic stress) acts as an additional 
source term for the gravitational potential IlSHH], changing the CMB anisotropy power 
spectrum at the level of ~ 20%. 

The delay of the epoch of matter-radiation equality shifts the LSS matter power 
spectrum turnover position toward larger angular scales, suppressing the power at 
small scales. In particular, the non-zero neutrino chemical potential leads to changes 
in neutrino free-streaming length and neutrino Jeans mass due to the increase of the 
neutrino velocity dispersion [IS, HE] . 

After WMAP3 data release there are many works aiming to constrain N^ff from 
cosmological observations [HI [Ml [371 113 HHl [12]. Their results suggest large values 
for Ngff within 95% CL interval, some of them not including the standard value 3.046 
[m [311 EZ] • Recently Ref . [50] argues that the discrepancies are due to the treatment of 
the scale-dependent biasing in the galaxy power spectrum inferred from the main galaxy 
sample of the Sloan Digital Sky Survey data release 2 (SDSS-DR2) [511 [52] and the large 
fluctuation amplitude reconstructed from the Lyman— a forest data [53] relative to that 
inferred from WMAP 3. 

Discrepancies between BEN and cosmological data results on iVe// was interpreted as 
evidence of the fact that further relativistic species are produced by particles decay 
between BEN and structure formation [48l[49]. Other theoretical scenarios include the 
violation of the spin-statistics in the neutrino sector [51], the possibility of an extra 
interaction between the dark energy and radiation or dark matter, the existence of a 
Erans-Dicke field which could mimic the effect of adding extra relativistic energy density 
between EEN and structure formation epochs [55] . 

A lower limit to iVe// > 2.3 (95% CL) was recently obtained from the analysis of 
the WMAP 5-year (WMAP5) data alone [56j, while the combination of the WMAP5 
data with distance information from baryonic acoustic oscillations (EAO), supernovae 
(SN) and Hubble constant measured by Hubble Space Telescope (HST), leaded to 
Neff = A A ±1.5 (68% CL), fully consistent with the standard value [57] . 



The extra energy density can be splitted in two distinct uncorrelated contributions, 
first due to net lepton asymmetry of the neutrino background and second due to the 
extra contributions from other unknown processes: 

^Kff = ^Kff{i)+^N:}). (2) 

The aim of this paper is to obtain bounds on the neutrino lepton asymmetry and 
on the extra radiation energy density by using WMAP5 data in combination with 
most of the existing CME and LSS measurements and self-consistent EEN priors on 
Yp. We also compute the sensitivity of the future Planck experiment [58] for these 
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parameters testing the restrictions on cosmological models with extra relativistic degrees 
of freedom expected from high precision CMB temperature and polarization anisotropy 
measurements. 

2. Leptonic asymmetric cosmological models 

The density perturbations in leptonic asymmetric cosmological models have been 
discussed in literature [591 EOl HSl HHl [631 [M] • We applied them to modify the Boltzmann 
Code for Anisotropies in the Microwave Background (CAMB) [65l|66l|67j to compute the 
CMB temperature and polarization anisotropies power spectra and LSS matter density 
fluctuations power spectra for the case of three degenerated neutrinos/antineutrinos with 
the total mass m^ and degeneracy parameter ^^. As neutrinos reach their approximate 
chemical potential equilibrium before BBN epoch [311 [321 [22], we consider in our 
computation that all three flavors of neutrinos/antineutrinos have the same degeneracy 
parameter ^^. 

When the Universe was hot enough, neutrinos and antineutrinos of each flavor z/j behave 
like relativistic particles with Fermi-Dirac phase space distributions: 

^^M = ^E^jn-i^ ^ ^ ^ ^M = ^E-.jT.-i, ^ 1 (^ = e,/i,r) (3) 



where Ey^ = ^J cp- + a^niy. is one flavor neutrino/antineutrino energy and g = ap is the 
comoving momentum. Hereafter, a is the cosmological scale factor (oq = 1 today). The 
mean energy density and pressure of one flavor of massive degenerated neutrinos and 
antineutrinos can be written as: 

P., + Pp. = (kBTy)' j^ ^ q^EyXJ'uM) + ^vM)) , (4) 

KPu. + PpJ = {ksT^Y r ^ ^{:FyM) + ^^M)) ■ (5) 

We modify in CAMB the expressions for the energy density and the pressure in 
the relativistic and non-relativistic limits for the degenerate case [16] and follow the 
standard procedure to compute the perturbed quantities by expanding the phase space 
distribution function of neutrinos and antineutrinos into homogeneous and perturbed 
inhomogeneous components [671 EHl [69]. Since the gravitational source term in the 
Boltzmann equation is proportional to the logarithmic derivative of the neutrino 
distribution function with respect to comoving momentum, dln(J^,^. + jFp.)/dlng, we 
also modify this term to account for ^jy 7^ [lEl 163] . 

As mentioned before, the BBN theory gives strong constraints on A^e// and ^^ by 
comparing the measured light element abundance with the theoretical predictions. The 
only free parameter is the baryon to photon ratio, tjb = ni,/ny, that is obtained from 
the determination of fi^/i^ from CMB measurements. 

In particular, the '^He mass fraction, Yp, affects the CMB angular power spectra through 
its impact on different evolution phases of the ionization/recombination history [70] . 
As previously demonstrated [711 221 [73], the impact of self-consistent BBN prior on 
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Yp has a net impact on parameter inference, improving the bounds on cosmological 
parameters compared to the analysis which treat 1^ as a constant or free parameter. 
We use the pubhc BBN code PArthENoPE [71] to compute the dependence of Yp on 
flbh"^, ANeff as given in equation (2) and $,^. The accuracy on Yp obtained by using 
PArthENoPE code is ~ 10~^, being only limited by the experimental uncertainty on 
the neutron lifetime [76]. We also modify the recombination routine Recfast vl.4 [75] 
of the CAMB code to account for Yp dependence on Q^h"^, AN^fj and C,^. 
In our computation we implicitly assume that the value of Ngff is not changed between 
the epoch of BBN and matter-radiation equality. 

3. Analysis 

We use the CoSMoMC Monte Carlo Markov Chain (MCMC) public package [77] 
modified for our extended 6+3 parameter space to sample from the posterior distribution 
giving the following experimental datasets: 

• The WMAP5 temperature and polarization CMB measurements [561 EZl [78] com- 
plemented with the CMB measurements from Boomerang [79l [80] , ACBAR [81] 
and CBI [82] experiments. 

• The LSS power spectrum of the matter density fluctuations inferred from the 
galaxy clustering data of the Sloan digital Sky Survey (SDSS) [SH [521 [83l [H] 
and Two-degree Field Galaxy Redshift Survey (2dFGRS) [HSj. In particular, the 
luminous red galaxies (LRG) sample from SDSS data release 5 (SDSS-DR5) [831 [M] 
has more statistical signiflcance than the spectrum retrieved from the SDSS main 
galaxy sample from data release 2 (SDSS-DR2) [5ll |52] eliminating the existing 
tension between the power spectra from SDSS-DR2 and 2dFGRS. For this reason 
we include in our analysis the matter power spectra from SDSS-LRG and 2dFGRS. 
We consider SDSS-LRG data up to kmax — 0.2h Mpc~^ and the 2dFGRS data up 
to kmax — 0.14h Mpc~^ and apply the corrections due to the non-linearity behavior 
and scale dependent bias [83|, connecting the hnear matter power spectrum, Pun^k), 
and the galaxy power spectrum, Pgai{k), through: 

P,al{k)=b^^-^^^Pun{k), (6) 

where the free parameters b and Qni are marginalized. 

• The luminosity distance measurements of distant Type la Supernovae (SNIa) 
obtained by Supernova Legacy Survey (SNLS) [86] and the Hubble Space Telescope 

m- 

• The BBN constraints on Yp as obtained from PArthENoPE code, allowing Qfyh"^ 
ANeff and ^^ to span the ranges indicated in Table 1. 
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Table 1. The free parameters of our model, their fiducial values used to generate the 
PLANCK-like simulated power spectra and the prior ranges adopted in the analysis. 



Parameter 


Symbol 


Fiducial value 


Prior range 


Baryon density 


n^h' 


0.022 


0.005 -^ 0.04 


Dark matter density 


Vtdmh 


0.11 


0.01 ^ 0.5 


Hubble constant 


Ho 


70 


40 ^ 100 


Redshift of reionization 


Zre 


11 


3^ 20 


Scalar spectral index 


ris 


0.96 


0.5 ^ 1.3 


Normalization 


ln[10^°v4; 


3.264 


2.7^4 


Neutrino density 


n,h^ 


0.01 


0^0.3 


Neutrino degeneracy parameter 


iu 





-1 -^ 1 


Number of extra rel. d.o.f. 


AN°f} 


0.046 


-3^3 


Helium mass fraction 


yp 


0.248 


0.07^0.6 



Hereafter, we will denote WMAP5+SDSS-DR5+2dFGRS+SNIa+BBN data set as 
WMAP5+A11. 

We perform our analysis in the framework of the extended ACDM cosmological 
model described by 6 + 3 free parameters: 

e = {n,h\ namh\ Ho, ^.e, n„ a, ^./^^ e., aat^J) . (7) 

y ^ / 

standard 

Here Qbh"^ and ^dmh"^ are the baryon and cold dark matter energy density parameters, 
Hq is the Hubble expansion rate, Zre is the redshift of reionization, n^ is the scalar 
spectral index of the primordial density perturbation power spectrum and Ag is its 
amplitude at the pivot scale /c* = 0.002 hMpc~^. The additional three parameters 
denote the neutrino energy density fij^/i^, the neutrino degeneracy parameter ^^ and 
the contribution of extra relativistic degrees of freedom from other unknown processes 
AN^lil . Table 1 presents the parameters of our model, their fiducial values used to 
generate the PLANCK-like simulated power spectra and the prior ranges adopted in the 
analysis. 

For the forecast from PLANCK-like simulated data we use the CMB temperature 
(TT), the polarization (EE) and their cross- corelat ion (TE) power spectra obtained 
for our fiducial cosmological model for multipoles up to / = 2000 and the expected 
experimental characteristics of the Planck frequency channels presented in Table 2 
[58] . We assume a sky coverage of fsky = 0.8. 

Following the method described in Ref . [6T| [62] , for each frequency channel we consider 
an homogeneous detector noise with the power spectrum given by: 

A^f, = (^feA„)2exp'('+i)^^'/«i°2 ce{T,P) (8) 
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Table 2. The expected experimental characteristics for the Planck frequency 
channels considered in the paper. At and Ap are the sensitivities per pixel for 
temperature and polarization maps. We assume a sky coverage of fsky = 0.8 



V 


FWHM 


At 


Ap 


(GHz) 


(arc-minutes) 


(/iK) 


(/xK) 


100 


9.5 


6.8 


10.9 


143 


7.1 


6.0 


11.4 


217 


5.0 


13.1 


26.7 



where v is the frequency of the channel, Q^, is the FWHM of the beam and Ac are the 
corresponding sensitivities per pixel of temperature (T) and polarization (P) maps. The 
global noise of the experiment is obtained as: 



t^t 



Ew 



Lu) 



(9) 



We assume uniform prior probability on parameters G (i.e. we assume that all values 
of parameters are equally probable) and compute the cumulative distribution function 
C{0) = fa C(Q)dQ/ L ""'"' CdQ , quoting as upper and lower intervals at 68% CL the 
values at which C{9) is 0.84 and 0.16 respectively. For the case of neutrino mass, m^, we 
quote the upper limit at 95%CL to facihtate the comparison with other measurements. 

4. Results 

We start by making a consistency check, verifying that by using WMAP5+A11 data and 
imposing ^u = ^ and AN°'jh = priors we obtain results in agreement with the ones 
obtained by WMAP collaboration [56|[57]. 

In order to understand how the extra relativistic energy density and the leptonic 
asymmetry affect the determination of other cosmological parameters, we compute 
first the likelihood functions for WMAP5+A11 by imposing C,u = ^ prior, extending 
then our computation over the whole parameter space for WMAP5-I-A11 and Planck- 
like simulated data. In Figure 1 we compare the marginalized likelihood probabilities 
obtained for the main cosmological parameters. 

As neutrinos with eV mass decouple when they are still relativistic {Tdec ~ 2 MeV), 
the main effect of including AA^g// is the change of relativistic energy density. This 
changes the redshift of matter-radiation equality, z^q, that affects the determination of 
flmh'^ from CMB measurements because of its linear dependence on N^fj [57] : 

Qruh"^ 1 

^^"^'1= n,h? 1 + 0.2271Ar,^/ • ^^°^ 
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Figure 1. The marginalized posterior likelihood probabilities of the main cosmological 
parameters obtained for: WMAP5+AU with AN°j'j = and ^^ = priors (black hnes), 
WMAP5+A11 with ^j, = prior (green hnes), WMAP5+A11 with no neutrino priors 
(blue lines) and PLANCK-like simulated data with no neutrino priors (red lines). 



Here Q^h'^=2A69 xlO^^ is the present photon energy density parameter for 
Tcmb = 2.725 K. As a consequence, N^jf and Vt^nh? are hnearly correlated, with the 
width of degeneracy hne given by the uncertainty in the determination of Zeq- 
In Figure 2 we compare the marginahzed posterior hkehhood probabihties of neutrino 
parameters, Yp and Zf,q obtained in our models. The mean values of these parameters 
and the corresponding (68% CL) error bars are given in Table 3. 

The LSS measurements provides an independent constraint on Qmh"^ which helps to 
reduce the degeneracy between this parameter and Zeq. From WMAP5+A11 data with 
AN°j'j = and ^;^ = priors we find a mean value of z^q = 3158 ± 68 (68% CL). One 
should note that mean value of Zeq for the standard ACDM model with N^ff = 3.046 is 
z,q = 3176+}^^ (68%CL) from WMAP5 data only [57|. 

Figure 3 presents the joint two-dimensional marginalized distributions (68% and 



95% CL) in fimh"^ - Nejf plane (left panel) and Vtmh^ 



^eq 



plane (right panel). The 



thick solid lines in the left panel show the 68% and 95% CL limits calculated from the 
corresponding limits on Zeq obtained form the WMAP5+A11 with S,v = ^ prior by using 
equation (10). 

When we transform A^e// axis of the left panel to z^q axis from right panel we observe a 
strong degeneracy between Zeq and VLmh?- This is valid for both WMAP5+A11 with ^j,=0 
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Figure 2. The marginalized posterior likelihood probabilities of the neutrino mass, 
TO,y, degeneracy parameter, ^^, leptonic asymmetry, Lj^, helium mass fraction, Yp, 
and redshift of matter-radiation equality, z^q, for: WMAP5-I-A11 with AN°j^r — and 
£,^ = priors (black lines), WMAP5-I-AU with ^^ = prior (green lines), WMAP5+A11 
with no neutrino priors (blue lines) and PLANCK-like simulated data with no neutrino 
priors (red lines). 



prior and WMAPS+AU with no neutrino priors. For the last case, the non-zero neutrino 
chemical potential augments the value of Neff by ANeff{^u) as given in equation (1). 
This imply a larger expansion rate of the Universe, an earlier weak process freeze out 
with a higher value for the neutron to proton density ratio, and thus a larger value of 
Yp. On the other hand, a non-zero value of the electron neutrino chemical potential 
shifts the neutron-proton beta equilibrium, leading to a large variation of Yp [721 [73] . 
Left panel from Figure 4 presents the two-dimensional marginalized joint probability 
distributions (68% and 95% CL), showing the degeneracy between Yp and A^e//- The 
total effect of ^,y 7^ is a noticeably increase of the projected error on Neff. 

As the anisotropic stress of neutrinos leaves distinct signatures in the CMB power 
spectrum which are not degenerated with Qmh'^, we conclude from this analysis that we 
observe a non-zero value of A^e// from WMAP5-I-A11 data mainly via the change of z^q 
rather than by the effect of neutrino anisotropic stress. 

Our conclusion is not in agreement with the claim of WMAP team concerning 
the strong evidence of neutrino anisotropic stress from a similar analysis of the 
WMAP5+BA0+SN+HST data [57]. 

We obtain from our analysis a mean value of A'e// = 3.026^Q;ggo (68% CL) from 
WMAP5+A11 with ^^ = prior and N^ff = 3.256toiZ (68% CL) form WMAP5+A11 
with no neutrino priors. These values bring an important improvement over the similar 
result obtained from WMAP5+BA0+SN+HST data: N^ff = 4.4 ± L5 (68% CL). 
From PLANCK-like simulated data the 68% error is a{Neff) ~ 0.3. 

The analysis of WMAP5-I-A11 data with iu = prior provides a strong bound on 
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Figure 3. Two-dimensional marginalized joint probability distributions (68% and 
95% CL) showing the degeneracy between: ftmh^ and N^ff (left panel) and flmh^ and 
Zeq (right panel). The thick solid lines in the left panel show the 68% and 95% CL 
limits calculated from the corresponding limits on z^q obtained form the WMAP5+A11 
with ^^ = prior by using equation (10). The contours show: WMAP5-I-A11 with 
AN°j'j = and ^t, = priors (black hues), WMAP5+AU with ^^ = prior (green 
lines), WMAP5-I-A11 with no neutrino priors (blue lines) and PLANCK-like simulated 
data with no neutrino priors (red lines). 



helium mass fraction oi Yp = 0.2486 ± 0.0085 (68% CL), that rivals the bound on Yp 
obtained from the conservative analysis of the present data on helium abundance [2]. 
Under the assumption of degenerated BBN this bound is weakened, leading to 



Yr. 



0.2487 



+0.0451 
-0.0484 



(68% CL) from WMAP5-I-A11 with no neutrino priors, that reflects 
the strong dependence of Yp on ^i,. 

From PLANCK-like simulated data the 68% error on Yp is o"(l^) = 0.0133, fully consistent 
with Yp bounds obtained from the conservative analysis of the present data on helium 
abundance. 

We get for neutrino degeneracy parameter a bound of —0.216 < ^jy < 0.226 
(68%CL) that represent an important improvement over the similar result obtained 
by using the WMAP 3-year data [i5] . 

The CMB only is able to constrain C,u through its contribution to the radiation energy 
density during radiation domination epoch and the BBN constraints on Yp. 
We find that the sensitivity of CMB to Yp from PLANCK-like simulated data will reduce 
the error on C,u down to cr{^u) — 0.089 (68% CL), value fully consistent with the BBN 
bounds. 

In the right panel from Figure 4 we compare the two-dimensional marginalized joint 
probability distributions (68% and 95% CL) between iVe// and ^^, showing the 
potentiality of the future high sensitivity CMB measurements to reduce the degeneracy 
between these parameters. 

The inclusion in the analysis of LSS data reduces significantly the upper limit of 
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Figure 4. Two-dimensional marginalized joint probability distributions (68% and 
95% CL) showing the degeneracy between: Yp and Ngff (left panel) and N^/f and 
^i, (right panel). The contours show: WMAPS+AU with S^i, — prior (green lines), 
WMAP5+A11 with no neutrino priors (blue lines) and PLANCK-like simulated data 
with no neutrino priors (red lines). 



the neutrino mass to m^ < 0.284 eV (95% CL) from WMAP5+A11 with AN^jf = and 
^^ = priors. One should note that the analysis of WMAP5+BA0+SN+HST leaded 
to m^ < 0.61 eV (95% CL) [57]. 

Under the assumption of degenerated BBN this bound is weakened due to the degeneracy 
between A^e// and C,u- From WMAP5-I-A11 with no neutrino priors we find rui, < 0.535 
eV (95%CL), while from PLANCK-like simulated data with no neutrino priors we obtain 
m^ < 0.440 eV. 



5. Summary and conclusions 

In this paper, we set bounds on the radiation content of the Universe and neutrino 
properties by using the WMAP-5 year CMB measurements complemented with most of 
the existing CMB and LSS measurements (WMAP5-I-A11), imposing also self-consistent 
BBN constraints on the primordial helium abundance, which proved to be important in 
the estimation of cosmological parameters [72l [73] . 

We consider lepton asymmetric cosmological models parametrized by the neutrino 
degeneracy parameter C,u and the variation of the relativistic degrees of freedom, AN°'jh, 
due to possible other physical processes occurred between BBN and structure formation 
epochs. 

From our analysis we get a mean value of the effective number of relativistic 
neutrino species form WMAP5-I-A11 with no neutrino priors, of N^jf = 3.256lo'g4^ 
(68% CL), bringing an important improvement over the similar result obtained from 
WMAP5+BA0-FSN+HST data [57]. 
Although the LSS measurements provides an independent constraint on Qmh'^ which 
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helps to reduce the degeneracy between this parameter and the redshift of matter 
radiation equahty, Zeq, we find a strong correlation between Q^h"^ and z^q, showing that 
we observe a non-zero Ngfj value from WMAP5+A11 data mainly due to the change of 
Zeq, rather than by the neutrino anisotropic stress. 

The analysis of WMAP5+A11 data with C,u = ^ prior provides a strong bound on 
helium mass fraction of 1^ = 0.2486 ± 0.0085 (68% CL), that rivals the bound on Yp 
obtained from the conservative analysis of the present data on helium abundance [2]. 
Under the assumption of degenerated BBN this bound is weakened to Yp = 0.2487lo;o4g4 
(68% CL), refiecting the strong dependence of Yp on ^^. 

From the analysis of WMAP5+A11 data we get for neutrino degeneracy parameter 
a bound of -0.216 < ^^ < 0.226 (68%CL). 

The inclusion in the analysis of LSS data reduces the upper limit of the neutrino 
mass to nil, < 0.535 eV (95% CL), from WMAP5+A11 with no neutrino priors, with 
respect to the values obtained from the analysis from WMAP5-only data [56] and 
WMAP5-I-BA0+SN+HST data [57]. 

The analysis of WMAP5+A11 measurements bring also an important improvement 
over the similar results obtained by using WMAP 1-year measurements complemented 
with LSS data [60] and the WMAP 3-year data alone [l5l[88]. 

We forecast that the CMB temperature and polarization measurements observed 
with high angular resolutions and sensitivity by the future Planck satellite will reduces 
the errors on ^^ and Yp down to a{^^) ~ 0.089 (68% CL) and a{Yp) = 0.0133 (68% CL) 
respectively, values fully consistent with the BBN bounds on these parameters [T2] . 
Our forecasted errors on the cosmological parameters from PLANCK-like simulated data 
are also consistent those obtained in Ref. 
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